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Abstract 

The long-standing problem of the effect of correlations on the fer- 
romagnetism of transition metals is apparently nearing solution. The 
ferromagnetism of transition metal compounds, for instance doped 
manganites, poses a new question: is there some kind of orbital order 
coexisting with itinerant ferromagnetism? The ideas and techniques 
introduced by Gutz wilier should be of use again. 

1 Introduction 

In the years 1963 to 1965, Gutzwiller published three highly influential papers 
[H, ^ ^ about the "Effect of Correlation on the Ferromagnetism of Transition 
Metals" . The problem had been long considered before but it came to a dead- 
lock: mean fleld criteria (Stoner theory) gave ferromagnetism all too easily. 
Since Van Vleck's classic analysis it was also known what the difficulty 
is: Stoner theory neglects the correlation between opposite-spin electrons, 
and therefore overestimates the kinetic energy gain which arises from mak- 
ing the spins parallel. One needed a method for describing local correlations 
between itinerant electrons. Clearly, the problem is ubiquitious in electronic 
systems, and has many manifestations beside the onset of ferromagnetism. 
If one can say why Fe, Co, and Ni are ferromagnetic while the corresponding 
Ad elements Ru, Rh, and Pd are not, one can probably also say why CoO 
is insulating - perhaps even more easily. And so it was to turn out. The 
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threefold contribution of Gutzwiller: a hamiltonian, a variational method, 
and an approximate analytical solution, proved to be of universal value; and 
of all the related problems, the ferromagnetism of the transition metals was 
one of the most difficult. 

In the next Section, I briefly describe why the problem of ferromagnetism 
is particularly difficult, and what our present view is. Though it does not 
come into this paper, it may be mentioned that the puzzle of the high-Tc 
superconductivity of cuprates inspired a keen re-examination of all correla- 
tion problems, and thus indirectly contributed to a renaissance of the subject 
pioneered by Gutzwiller. Lately, we have been witnessing the important de- 
velopment that the combination of the Gutzwiller variational approach with 
LDA band structure calculation promises to become one of the major tech- 
niques of computational solid state physics. It has been demonstrated that 
starting with nine (s, p, d) orbitals per site, and using a judiciously chosen 
number of variational parameters to describe the possible on-site spin-spin, 
orbital-spin and orbital-orbital correlations, an impressive array of the elec- 
trical, magnetic and optical properties of nickel can be calculated, with an 
overall quality which surpasses that of spin density functional theory [0]. 
One might perhaps claim that the program initiated by Gutzwiller has been 
brought to a successful conclusion, using the very tools he devised. The at- 
tractive feature of the approach outlined in 0| would be that (in spite of 
the inevitable complexity of the procedure) one still has a clear idea which 
correlations have been taken into account and what their role is. 

Even within the family of "id systems, the nature of ferromagnets is quite 
varied and while essentially full undestanding may be in sight for some of 
them, basic questions remain to be answered for others. The ferromagnetic 
transition metal compound Lai-^^Sr^^MnOs comes to mind f^. One may as- 
sume that if we understand the ferromagnetism of Ni, we should understand 
that of Lai_a;Sra,Mn03 even better, since this compound gives a better real- 
ization of a Gutzwiller-Kanamori-Hubbard model with partial band filling 
than Ni does: it does not have the complication of an overlapping 4s band. 
However, the nature of ferromagnetic manganites is a matter of lively debate 



[|, |T0|, |TT|. Though the problem of complex orbital order (Sec. 3) does not 
have the same broad significance as ordinary ferro- or antiferromagnetism, it 
has some amusing parallels with the long-studied problem of spin ordering, 
and leads to a series of questions which are not unlike those which Gutzwiller 
set out to investigate long ago. 
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2 Ferromagnetism of Correlated Electrons 



In [U and 0, Gutzwiller discussed the key ingredients which give rise to 
the ferromagnetism of transition metals, particularly nickel. He distingushed 
between conduction (4s), and the almost localized 'valence' {3d) electrons, 
whereby the interactions among the latter kind of electrons were identified 
as the cause of magnetism. Though his description in words amounts to a 
recipe for the periodic Anderson model (including hybridization), Gutzwiller 
decided to omit the conduction electrons, and kept only the band of inter- 
acting electrons in the HamiltonianQ. First, let us write down the 1-band 
version 

Hi_band = -tJ2 E(4cja + H.c.) + "-JT^ji • (1) 

(iJ) j 

The dichotomy of itinerant and localized features (embodied in the hopping 
and on-site interaction terms, resp.) makes the analysis of this simple-looking 
hamiltonian extremely difficult. Being aware of the pitfalls of a mean- field- 
like smearing-out of the interaction term, Gutzwiller insisted that the many- 
electron state should contain local correlations. This is achieved by gradually 
suppressing the double occupation of sites in a spin-polarized Fermi sea: 

|fc|<fcF-|- |fc|<fcp| 

i^)=n[i-(i-^)^jT^ji] n 4t n 4jo) (2) 

j k k 

There are two variational parameters: rj which controls the local correlations, 
and the polarization density m = (N^—Ni) / L = n^—rii, where L is the num- 
ber of lattice sites. ^ can be used for any band filling n = + rti. In spite 
of its apparent simplicity, ^ is very difficult to handle: in order to evalute 
rid, the density of doubly occupied sites, the Fermi sea has to be decomposed 
into ~ localized configurations. The full implementation of the varia- 
tional procedure became possible only 24 years later [|12|. However, a great 
variety of approximate variational solutions came earlier, thanks to a tech- 
nical innovation announced in Gutzwiller's third paper 0: the Gutzwiller 
approximation]^ (in short GA). It proposes a practicable analytical method 
for solving the many-electron variational problem formulated in the letter 

^Nowadays we refer to it as the Hubbard model, but it was introduced independently 
by Gutzwiller Hubbard and Kanamori in the same year. 

^Ref. IQ gives an introduction to the Gutzwiller variational method. 
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[|^]. The essential result is that the kinetic energies get scaled down from 
their U = values by the spin- dependent renormalization factors qa, and the 
energy expression becomes 



(^l^l-bandl^) 



= l^^'^ + Und (3) 

GA o- 



Here q„ still depends on which has taken over the role of the original vari- 
ational parameter rj. It is enlightening to quote the strong coupling limiting 
form of 

q^{U ^ oo) = ^—^ . (4) 

For a fully polarized system = n, = 1, the kinetic energy keeps its 
bare value. This echoes the old exchange hole argument for the arising of 
ferromagnetism: parallel-spin electrons can avoid each other, there is no 
interaction energy penalty, but it comes at a price: the majority spin band 
has to be filled up to a higher Fermi energy. In contrast, the paramagnetic 
solution keeps the Fermi level deeper down in the band, but the kinetic energy 
gain is reduced by correlation-induced band narrowing: 

1 — 11 

q^{U ^ oo) = q^{U ^ oo) = ^—^^ < 1 . (5) 

For bands with a symmetrical density of states, the ferromagnetic and param- 
agnetic solutions reach a balance at about quarter filling, and the low-density 
system remains non- magnetic even at t/ ^ oo [Q. This should be contrasted 



with the Stoner result that the system becomes ferromagnetic at arbitrary 
filling if U is large enough. The criteria for the onset of antiferromagnetism 
are also strongly modified [0, 0]. 

In fact, at half-filling {n = 1), q{U oo) oc (1—"^) = 0, signalling that we 
are within the Mott insulating phase. Being primarily interested in itinerant 
ferromagnetism, Gutzwiller did not pay attention to the particular case of 
integral band filling but it was soon realized by Brinkman and Rice []TB| that 



a critical divergence of the effective mass of the electrons on the metallic 
side of the Mott transition is a corollary of the results presented in The 
increasing "heaviness" of the fermions is manifested in an enhancement of the 
electronic specific heat and the spin susceptibility, and a reasonable Wilson 
ratio can be derived. This was the beginning of our understanding of the 
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T = Mott transition as a quantum phase transition; up-to-date studies 
using dynamical mean field theory still yield results which bear remarkable 
similarity to the Gutzwiller-Brinkman-Rice scenario (see, e.g. ||17|). 

Subsequent studies of the periodic Anderson model (in which the strongly 
correlated band has non-integral filling) suggested a similar interpretation of 



the heavy Fermi sea of almost-integral- valent /-electron systems [|18| , |T^]. 
For a wide range of problems, the Gutzwiller method allowed to develop a 
pictorial way of thinking about local correlation effects, and also provided 
the technique to do the corresponding calculations. The generalization of 
the ideas to degenerate bands, and to several correlated bands, takes a lot of 



work but is conceptually straightforward ||2T| , p2| . The use of the variational 
method is not constrained by the GA; we now understand that GA is really 
an infinite-dimension [D = oo) approximation, and we can go beyond it by 
either exact evaluation []T^, or by calculating 1/D-corrections for D > 2 [^]. 
Though the Gutzwiller method goes a long way to correct the Hartree- 



Fock results for the magnetic phase diagram, we have reason to suspect 
that it still overestimates the chances of ferromagnetic ordering. It is a typical 
Gutzwiller result that the extent of the FM phase along the n axis is largest 
at U ^ oo. This justifies another look at the mechanism driving ferromag- 
netism in the strong coupling limit. Performing the canonical transformation 
which eliminates double occupation in lowest order, we arrive at the familiar 
t-J model: 

= -^EE(5U. + H.c.) + ^Efs.-s,-^) 

{id) ^ \ / 

+ (three — site terms) (6) 

Let us note that the 1-band model generates only one kind of effective spin- 
spin interaction and that is antiferromagnetic. This need not rule out the 
possibility of a ferromagnetic ground state but if it exists, the reason for it is 
not straightforward. Given that itinerant ferromagnetism (if any) does not 
originate from the interaction term in (|^), it has to arise from the projected ki- 
netic energy; it is an elusive, and often fragile, phenomenon. Since the nature 
of a uniformly spin-polarized state is so easy to grasp, it is somewhat surpris- 
ing to find out that the ferromagnetism of the 1-band model poses a much 
more difficult many-body problem than antiferromagnetism, and that even 
for three-dimensional systems, mean field predictions are grossly misleading 
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2^. After thirty years of continuing efforts we know that Gutzwiller's imme- 
diate reaction [|l| : "It seems rather that the exact ground state of the Hamil- 
tonian (1) is never ferromagnetic" was not completely right but it captured 
the essential point that the existence of a robust ferromagnetic phase is not a 
generic property of the 1-band model. For bipartite lattices with a symmet- 
rical density of states, even if a ferromagnetic phase exists, it is only at such 
large coupling strengths that we can safely forget about it in connection with 
real systems. However, intermediate-coupling ferromagnetism was shown to 
appear for non-bipartite lattices (e.g. fee) or/and a non-symmetric band den- 
sity of states |2^. The relevance of these factors was already pointed out in 
the pioneering works of Gutzwiller, Kanamori, and Hubbard; Gutzwiller also 
realized that direct exchange which is ordinarily negligible compared to other 
Coulomb terms, can become important when the magnetic and non-magnetic 
states are finely balanced. 



3 Orbital Liquid vs Complex Orbital Order 

Looking back, we might say that the detour via the single-orbital model (|^) 
made the understanding of ferromagnetism unnecessarily difficult. In his 
first papers, Gutzwiller clearly stated that the degeneracy of the d-band, and 
the intra-shell Hund exchange, are of vital importance. However, technical 
difficulties forced us to start with the 1-band model. This has brought many 
important results (foremost an understanding of the Mott transition), but 
it gradually became clear that the model is very reluctant to yield itinerant 
ferromagnetism at reasonable couplings. We now largely understand why 
this is the case, but perhaps we should never have worried too much: most 
of the itinerant ferromagnets we are interested in are based on degenerate 
shells. Models with orbital degeneracy have several efficient mechanisms to 
give spin ferromagnetism. However, they bring the additional complications 
of orbital ordering or/and orbital fluctuations. 

A model with twofold orbital degeneracy is relevant for the doped ferro- 
magnetic manganites such SrajMnOs which are Mn^+-Mn'^"'" valence 
fluctuators. Both Mn^+ and Mn^"*" have stable T| cores in the 5 = 3/2 
high-spin state. The metallic nature arises from the fourth electron of Mn^+ 
ions which may hop around in the ii^-shell states. For a realistic descrip- 
tion of manganites, it is important to remember the T2 cores, and the T2-E 
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intra-atomic exchange |1^], but the basic questions of spin magnetism 
and orbital order can be discussed within the 2-band modeL We write it 
down for the pair of sites 1, 2, which is sufficient for deriving the strong- 
couphng effective hamiltonian. a and h denote the orbitals. In general, both 
the interorbital and intraorbital hopping amplitudes differ from zero 

^hop = -ta Y^{ClaaC2aa + H.C.) - tfc Y^{c\i,„C2ha + H.C.) 

(T a 

-tab IZ(clb^C2aa + c\^„C2ba + H.C.) . (7) 
a 

We do not write down the familiar interaction terms which contain the 
intra- and interorbital Hubbard terms (with Ua-, Ub, and Uab), and the ferro- 
magnetic Hund coupling —2 J. Actually taking the cubic E doublet would 
require a specific relationship between the parameters, but we keep the nota- 
tion general because we intend to apply it to the discussion of other crystal 
field levels as well. 

The 2x2 (spin x orbital) degeneracy of the shell states is conventionally 
expressed by representing the orbitals as r = 1/2 pseudospin states: ipa —>■ 

= +(1/2)) and (pb = —(1/2))- This allows to write the large-t/ 

effective hamiltonian as mixed spin-pseudospin interaction. Its derivation 
from the 2-band Hubbard model is straightforward but tedious. We do not 
quote the full result (the familiar Kugel-Khomskii model [^), merely the 
term which acts on spin triplets: 

^eff = |^-^(^ + 81-82) {constant + 4(t,t,- 4) rrT2 

+ 2ita - tbfrlT^ + Kb^ir^ + 4U(ta - 4)(rM + r^^^)} (8) 

The spin- singlet term is similar. Let us observe that in general, the interac- 
tion is anisotropic in pseudospin space. 

The usual choice of basis for the cubic E doublet is shown in Fig. The 
real basis states have lobes in specific directions. The hopping amplitudes in 
(|^ are not merely orbital index dependent, but depend also on the direction 
of the hopping. 

Let us now look for the ferromagnetic ground state of the two-band model. 
For the sake of simplicity, we may assume full spin polarization, so that we 
once again have only two states per site (the two orbitals), like in the simple 
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Figure 1: Left and middle: the conventional real basis functions ipa ~ — and 
ifb ~ — for the cubic E doublet. Right: the complex combinations fa i iifb 
are octupolar eigenstates. 

1-band model (0). It may seem that the remaining problem of an orbitally 
correlated state is at the same level of difficulty as the thoroughly stud- 



ied correlation problem of the ordinary Hubbard model However, the 

orbital index dependence of the hopping matrix elements brings additional 
complications: it is as if in the single-orbital model, we had spin-dependent 
hopping, with spatial and spin anisotropy. In the case of (|l|), if we find a 
spin-ferromagnetic state, we know that the total spin can point in any direc- 
tion. For the two-band model, polarization in the pseudospin ^r-direction is 
something totally different from polarization along the y-direction. 

The ground state we seek should be analogous to (Q) but how to specify 
the orbitals? Postulating an orbital order of the kind found in insulating 
LaMnOa would mean that each unit cell has an asymmetric form of the 
3d electron cloud, and an accompanying Jahn-Teller distortion. However, 
experiments tell us that in the ferromagnetic metal the Jahn-Teller distortion 
is absent^ and the unit cells should be thought of as cubic. An obvious idea 
is that this is the result of time-averaging: the ground state is orbitally 
disordered ("pseudospin paramagnetic"). Let us take the non-interacting 
two-band Fermi sea and Gutzwiller-project it (now two electrons can share 



^At least in some of the manganites. 
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the same site only if they are in different orbitals) 



\k\<kpa \k\<kpg 

i$)=n[i-(i-^)^jaT^M] n n ^/.tIq) (q) 

j k k 

where a, (3 refer to the two bands (the band structure is shown in |T^). In 
spite of orbital disorder, the ground state is unique (an orbital liquid ||TT|), 
in the same sense that the ordinary Fermi sea is a spin liquid. 

Is there some other possibility? Can we prescribe some kind of orbital or- 
der and still get a similar picture of a cubic ferromagnetic metal? Obviously, 
it could not be done by ascribing one of the real wave functions — y^) or 
\?)Z^—r'^) to the sites. However, there is another possibility. The order param- 
eter space is described by the product representation E®E = Ai + A2 + E. 
E belongs to the familiar quadrupolar doublet of (x^ — ?/^)-like and (3z^ — r^)- 
like order parameters, but we have not yet considered the single component 
order parameter with A2 symmetry. Since it transforms like xyz, it must be 
the magnetic octupolar moment L^LyL^ ||10| where the bar indicates sym- 
metrization over six terms 

L.LyL, ^ ^{L.LyL, + L,L,Ly + ...)= '-{L+L'L+ - L-UL') cc . (10) 

The = ±(1/2) eigenstates are the complex combinations 

\C±) = ^{\x'-y')±r\^z'-r^)) (11) 

First of all, let us observe that the charge density contours of these com- 
plex orbitals show cubic symmetry (Fig. |l|, right). Postulating the uniform 
octupolar state 

l^o) = n(iT+^6T)|0) (12) 
j 

it would show the cubic symmetry demanded from the ferromagnet. Nat- 
urally, would be a fully localized state with one electron at each site, 
so it is not yet what we want. To have a uniformly octupolar itinerant spin 
ferromagnet, we have to construct a trial state analogous to (^, a Gutzwiller- 
projected Fermi sea with a finite octupolar polarization . 



The arising of complex order as a possibility distinct from real order fol- 
lows from the lack of rotational symmetry in isospin space. There is another 
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interesting aspect: the possibility of time reversal invariance breaking with- 
out (the usual kind of) magnetism. Actually, it occurs whenever we have 
an orbital (non-Kramers) doublet to start with. The local basis can always 
be chosen in a complex form, and there will be a corresponding complex 
operator to which these are eigenstates (in the present case, it is the purely 
imaginary r^). But a complex quantity changes to its conjugate under time 
reversal, thus the state is breaking that symmetry. 

In which form time reversal invariance breaking happens, depends on 
the kind of doublet we start from. In the present case, it is broken by a 
magnetic octupolar component, which has been reintroduced into solid state 



physics only recently, mainly with reference to CeBg P5| , P6[ . It is, however, 
interesting to recall that the very problem we are discussing, namely the 
ordering patterns supported by cubic E subshells, has been considered in 
the early work by Korovin and Kudinov [^]. These authors remarked that 
postulating a state like "...leads to nonzero values of the third moments 
/ XaXi3X-yjs{r)dr of the current density j(r) in the crystal^ (the first moment, 
i.e., the usual orbital magnetic moment, is absent for the Eg representation)". 
The suggestion was then apparently forgotten because insulating LaMnOs 
has real orbital order (it is Jahn-Teller distorted), and the possibility of an 
essentially invisible orbital order in the cubic ferromagnetic phase was not 
taken seriously. However, the resurgence of general interest in higher mul- 
tipolar order, and the availability of sophisticated experimental techniques 
for its detection, led to a reconsideration of the octupolar phase |P, [T^. In 
particular, Takahashi and Shiba |10| determined the orbital analogue of the 
Penn phase diagram [^, comparing the threshold values for the RPA in- 



stabilities against various kinds of uniform and staggered orders, and found 
that staggered octupolar order is the first instability in a wide range of band 
filling. By its nature, theirs is a weak coupling calculation, and it is clearly 
desirable to complement it by a strong-coupling, Gutzwiller-like calculation. 

One reason why octupolar order may have long eluded consideration, is 
that it is not motivated by the strong coupling effective hamiltonian (P). For 
the cubic E doublet the hopping amplitudes in the x-y plane obey ta : h : 
\tab\ =3:1: v^S, thus the coefficient of rfrf vanishes. It is a little bit 



^The macroscopic magnetic octupolar moment introduced by Korovin and Kudinov is 
not literally the same as the order parameter ( p^ considered in , but the two concepts 
are related. 
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like the story of ferromagnetism in the 1-band model where we did not have 
ferromagnetic spin-spin interaction in the t-J model, but nonetheless felt 
justified in looking for a ferromagnetic ground state. 




Figure 2: Left and middle: real basis functions for the trigonal doublet E'tr- 
Right: the complex linear combinations have unquenched orbital moment along 
the trigonal z axis. 

It is interesting to compare the case of the usual E doublet with the 
doublet i?tr which is split off from the T2 level by a trigonal distortion, such 
as seen in the 3(i^ compound BaVSa ||3^ , |33| . The real wave functions are 
shown in Fig. ^ (left and middle), while the complex combination on the right. 
It apparently has a trigonal shape but can be shown to carry a remnant of 
unquenched orbital moment along the trigonal 2;-axis. The decomposition 
for the order parameter E®E = Ai + A2 + E is valid here, too, but the basis 
function of A2 is simply z ^ L^, thus time reversal invariance breaking 
can be realized with (orbital) magnetic dipole moment. 

It is a great honour to be permitted to dedicate this paper to Martin 
Gutzwiller on the occasion of his 75th birthday. 
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